OSCILLATION  OF  CYLINDERS  IN  OR  BELOW  THE  FREE  SURFACE  OF  DEEP  FLUIDS  Report  2375 


NAVAL  SHIP  RESEARCH  AND  DEVELOPMENT  CENT 


WASHINGTON,  D.C.  20007 


’>P  R  3  0.<L 


.y*> 

SO 

00 


flag 


OSCILLATIOi'l  OF  CYLINDERS  IN  OR  BELOW 
THE  FREE  SURFACE  OF  DEEP  FLUIDS 


i 


i 

i 


This  document  has  been  approved  ror  public  release  and  sale; 
its  distribution  is  unlimited. 


D 


MAY 


n 

0 


HYDROMECHANICS  LABORATORY 
RESEARCH  AND  DEVELOPMENT  REPORT 


October  1967 


Report  2375 


OSCILLATION  OF  CYLINDERS  IN  OR  BELOW  THE  FREE 
SURFACE  OF  DEEP  FLUIDS 


by 

W.  Frank 

ABSTRACT 

The  subject  of  cylinders  oscillating  in  or  below  the  free  surface  of  very 
deep  fluids  is  developed  as  a  boundary  value  problem  within  the  framework  of 
linear  free-surface  theory  by  distributing  source  singularities  over  the  sub¬ 
merged  portion  of  the  cylinders.  A  computer  program  for  calculating  the  hydro- 
dynamic  pressure,  force,  and  moment  on  these  cylinders  has  been  devised  by 
the  Naval  Ship  Research  and  Development  Center.  The  results  for  various  two- 
dimensional  shapes  are  given  in  graphical  form. 

ADMINISTRATIVE  INFORMATION 

This  study  was  authorized  by  the  Naval  Ship  Systems  Command  under  the  General 
Hydromechanics  Research  Program,  and  it  was  funded  under  Subproject  S-R009  01  01, 
Task  0100. 


TABLE  OF  CONTENTS 


Page 

ABSTRACT .  i 

ADMINISTRATIVE  INFORMATION .  i 

INTRODUCTION .  1 

FORMULATION  OF  PROBLEM .  2 

SOLUTION  OF  PROBLEM .  5 

LIMITATIONS  OF  GREEN’S  FUNCTION  -  TNTFJRAL 

EQUATION  METHOD .  9 

NUMERICAL  RESULTS  AND  CONCLUSIONS . 11 

ACKNOWLEDGMENTS .  18 

APPENDIX  A  -  EVALUATION  OF  PRINCIPAL  VALUE  INTEGRALS .  14 

APPENDIX  B  -  EVALUATION  OF  KERNAL  INTEGRALS .  18 

APPENDIX  C  -  EVALUATION  OF  POTENTIAL  INTEGRALS . 20 

APPENDIX  D  -  VELOCITY  POTENTIALS  FOB  VERY  SMALL  AND 

VERY  LARGE  FREQUENCIES .  24 

REFERENCES . .  89 


LIST  OF  FIGURES 


Figure  1  -  Polygonal  Approximation  to  Immersed  Part  of  Cylindrical 

Cross  Section  CQ .  26 

Figure  2  -  Kinematics  of  Roll  Oscillation .  26 

Figure  8  -  Adjoint  Interior  Boundary  Value  Problem  . .  26 

Figure  4  -  First  Irregular  Wave  Number  as  a  Function  of  B/T 

with  T  Fixed  for  Rectangular  Cylinders  . .  27 

Figure  5  -  Heave  Added  Mass  and  Damping  Coefficients  for 

Rectangular  Cylinder  . 27 

Figure  6  -  Scaled  Determinant  of  Matrix  of  Influence  Coefficients 

for  Heaving  Rectangular  Cylinder . 28 

Figure  7  —  Heave  Added  Mass  and  Damping  Coefficients  for 

Rectangular  Cylinder  .  28 

Figure  8  -  Heave  Added  Mass  and  Damping  Coefficients  for 

Semi-Immersed  Circular  Cylinder . 29 

Figure  9  -  Pressure  in  Phase  with  Displacement  on  Heaving 

Semi-Immersed  Circular  Cylinder  .  29 

Figure  10  —  Pressure  in  Phase  with  Velocity  on  Heaving 

Semi-Immersed  Circular  Cylinder  .  30 

Figure  11  —  Sway  Added  Mass  and  Damping  Coefficients  for 

Semi-Immersed  Circular  Cylinder .  30 

Figure  12  —  Pressure  in  Phase  with  Displacement  on  Swaying 

Semi-Immersed  Circular  Cylinder  . 31 

Figure  18  —  Pressure  in  Phase  with  Velocity  on  Swaying 

Semi-Immersed  Circular  Cylinder . 31 

Figure  14  -  Added  Mass  and  Damping  Coefficients  of  Fully  Submerged 

Heaving  Circular  Cylinder .  32 

Figure  15  -  Added  Mrss  and  Damping  Coefficients  of  Heaving 

Ogival  Cylinder .  32 

Figure  16  -  Pressur  e  in  Phase  with  Displacement  on  Heaving 

Ogival  Cylinder .  33 

Figure  17  —  Pressure  in  Phase  with  Velocity  on  Heaving 

Ogival  Cylinder .  33 

iii 


Page 

Figure  18  —  Heave  Added  Mass  and  Damping  Coefficients  for 

Midsection  of  Series  60,  Block  70  Ship . . .  34 

Figure  19  -  Sway  Added  Mass  and  Damping  Coefficients  for 

Midsection  of  Series  60,  Block  70  Ship . .  34 

Figure  20  —  Boll  Added  Moment  and  Damping  Coefficients  for 

Midsection  of  Series  60,  Block  70  Ship .  35 

Figure  21  -  Heave  Added  Mass  and  Damping  Coefficients  for 

Bow  Section  of  Series  60,  Block  70  Ship  .  35 

Figure  22  —  Sway  Added  Mass  and  Damping  Coefficients  for 

Bow  Section  of  Series  60,  Block  70  Ship  .  36 

Figure  28  -  Roll  Added  Moment  and  Damping  Coefficients  for 

Bow  Section  of  Series  60,  Block  70  Ship  . .  36 

Figure  24  —  Heave  Added  Mass  and  Damping  Coefficients  for 

Aft  Section  of  Series  60,  Block  70  Ship .  37 

Figure  25  -  Sway  Added  Mass  and  Damping  Coefficients  for 

Aft  Section  of  Series  60,  Block  70  Ship .  37 

Figure  26  -  Roll  Added  Moment  and  Damping  Coefficients  for 

Aft  Section  of  Series  60,  Block  70  Ship .  38 


iv 


NOTATION 


A**’*  Oscillation  amplitude  in  the  mth  mode 

B  Beam  of  cross  section  C0 

C0  Submerged  part  of  cross  sectional  contour  in  rest  position 

g  Acceleration  due  to  gravity 

/jf*'  Influence  coefficient  in  phase  with  displacement  on  the  ith 

midpoint  due  to  the  j  th  segment  in  the  m th  mode  of  motion 

Influence  coefficient  in  phase  with  velocity  on  the  »th 
midpoint  due  to  the  /  th  segment  in  the  nth  mode  of 
oscillation 

Mffi  Added  mass  force  or  moment  for  the  sith  mode  of 

oscillation  at  frequency  to 


N  Number  of  line  segments  defining  submerged  portion  of  half 

section  in  rest  position 

Damping  force  or  moment  for  the  si  th  mode  of  oscillation 
at  frequency  <o 

Direction  cosine  of  the  normal  velocity  at  t'th  midpoint 
for  mth  mode  of  oscillation 

P.V.  Cauchy  principal  value  of  integral 

pj")  Hydrodynamic  pressure  in  phase  with  displacement  for 

mth  mode  of  oscillation 

pj"*  Hydrodynamic  pressure  in  phase  with  velocity  for  mth 

mode  of  oscillation 

Q}m^  Source  strength  in  phase  with  displacement  along  /'th 

segment  for  mth  mode  of  oscillation 

C/7 Jv  Source  strength  in  phase  with  velocity  along  /  th  segment 

for  mth  mode  of  oscillation 

a  Length  variable  along  C0 

8-  /th  line  segment 

f  Draft  of  cross  section 
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Normal  velocity  component  at  tth  midpoint  for  mth  mode 
of  oscillation 

Abscissa  of  tth  midpoint 
Ordinate  of  tth  midpoint 
Ordinate  of  center  of  roll 

Complex  field  point  in  region  of  fluid  domain  (»  -  *  +  iy) 
Complex  midpoint  of  tth  segment  ( ai  -  +  iy.) 

Angle  between  f  th  segment  and  positive  x  axis 

Complex  variable  along  C0 

/'th  complex  input  point  along  C0 

Ordinate  of  /  th  input  point 

Wave  number  (  -  <u2/g) 

Ath  irregular  wave  number  for  adjoint  interior  problem . 
Abscissa  of  ;'th  input  point 
Density  of  fluid 

Velocity  potential  for  sith  mode  of  oscillation 

Radian  frequency  of  oscillation 

i.th  irregular  frequency  for  adjoint  interior  problem 
(Ath  eigenfrequency) 
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INTRODUCTION 


Hydrodynamic  research  of  horizontal  cylinders  oscillating  in  or  below  the  free  sur¬ 
face  of  a  deep  fluid  has  increased  in  importance  in  the  past  two  decades  and  has  been  stud¬ 
ied  by  a  number  of  investigators.  The  modern  history  of  this  subject  began  with  Ursell,1 
who  formulated  and  solved  the  boundary-value  problem  for  the  semi-immersed  heaving  circu¬ 
lar  cylinder  within  the  framework  of  linearized  free-surface  theory.  He  represented  the 
velocity  potential  *s  the  sum  of  an  infinite  set  of  multi  poles,  each  satisfying  the  linear  free- 
surface  condition  and  each  being  multiplied  by  a  coefficient  determined  by  requiring  the 
series  to  satisfy  the  kinematic  boundary  condition  at  a  number  of  points  on  the  cylinder. 

Grim2  used  a  variation  of  the  Ursell  method  to  solve  the  problem  for  two-parameter, 
Lewis-form  cylinders  by  conformal  mapping  onto  a  circle.  Tasai3  and  Porter,4  using  the 
Ursell  approach,  obtained  the  added  mass  and  damping  for  oscillating  contours  mappable 
onto  a  circle  by  the  more  general  Theodorsen  transformation.  At  the  Naval  Ship  Research 
and  Development  Center,  Ogilvie*  calculated  the  hydrodynamic  forces  on  completely  sub¬ 
merged  heaving  circular  cylinders. 

Despite  the  success  of  the  multipole  expansion-mapping  methods,  the  present  paper 
discusses  the  problem  from  a  different  view.  The  velocity  potential  is  represented  by  a 
distribution  of  sources  over  the  submerged  cross  section.  The  density  of  the  sources  is  an 
unknown  function  (of  position  along  the  contour)  to  be  determined  from  integral  equations 
found  by  applying  the  kinematic  boundary  condition  on  the  submerged  part  of  the  cylinder. 
The  hydrodynamic  pressures  are  obtained  from  the  velocity  potential  by  means  of  the  linear¬ 
ized  Bernoulli  equation.  Integration  of  these  pressures  over  the  immersed  portion  of  the 
cylinder  yields  the  hydrodynamic  forces  or  moments. 

A  simpler  approximation  to  toe  solution  of  tue  two-dimensional  hydrodynamic  problem 
was  used  in  toe  strip  theory  of  ship  motions  introduced  by  Korvin-Kroukovsky.6  The  solu¬ 
tions  of  two-dimensional,  water-wave  problems  for  ship  sections  by  multipole  expansion  and 
mapping  techniques  have  been  applied  to  this  strip  theory  by  Gerritama,7  Smith,8  Grim,9  and 
Vassilopoulos10  to  predict  the  heaving  and  pitching  motions  of  surface  vessels  moving  in 
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head  or  following  seas.  This  report  has  largely  been  motivated  by  the  desirability  of  devising 
a  computer  program,  based  on  strip  theory  and  independent  of  mapping  techniques,  to  predict 
the  response  of  surface  ships  moving  with  steady  forward  speed  in  oblique  as  well  as  bead 
or  following  seas  for  all  six  degrees  of  freedom. 

FORMULATION  OF  PROBLEM 

Consider  a  cylinder,  whose  cross  section  is  a  simply  connected  region,  which  is  fully 
or  partially  immersed  horizontally  in  a  previously  undisturbed  fluid  of  infinite  depth.  The 
body  is  forced  into  simple  harmonic  motion,  and  it  is  assumed  that  steady  state  conditions 
have  been  attained. 

The  two-dimensional  nature  of  the  problem  implies  three  degrees  of  freedom  of  motion. 
Therefore,  consider  the  following  three  types  of  osoillatory  motions:  vertical  or  heave,  hori- 
sontal  or  sway,  and  rotational  about  a  horizontal  axis  or  roll. 

To  use  linearised  free- surface  theory,  the  following  assumptions  are  made: 

1.  The  fluid  ie  incompressible  and  inviscid. 

S.  The  effects  of  surface  tension  are  negligible. 

8.  The  fluid  flow  is  irrotstional. 

4.  The  motion  amplitudes  end  velocities  are  ail  small  enough  that  all  but  the  linear 
terms  of  the  free  surface  condition,  the  kinematic  boundary  condition  on  the  cylinder,  and  the 
Bernoulli  equation  may  be  neglected. 

Far  complete  discussions  of  linearised  free-surfsoe  theory,  the  reader  is  referred  to  Stoker11 
and  Wehauaen  and  Laitone. 13 

Given  the  above  conditions  and  assumptions,  the  problem  reduces  to  the  following 
boundary-value  problem  of  potential  theory.  The  cylinder  is  forced  into  the  simple  harmonic 
motion  cos  (erf)  with  prescribed  radian  frequenoy  <u,  where  the  superscript  msy  take  on 
the  values  8,  8,  and  4  denoting  swaying,  heaving,  and  rolling  motions,  respectively.  It  is  re¬ 
quired  to  find  a  velocity  potential 

♦<«>  (,,  y.  t)  .  St  (*<">  <«,  y)  «-  <««)  (11 

satisfying  the  following  conditions: 

1.  The  Laplace  equation 


+  18) 

in  the  fluid  domain,  i.e.,  for  y  <  0  outside  the  cylinder; 
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2.  The  free  surface  condition 


*ttm)  +  t*jm)-0  w 

on  the  free  surface  y  -  0  outside  the  cylinder,  g  denoting  the  acceleration  of  gravity; 
8.  The  “bottom"  condition 


W*(m)\  -0  [4] 

4.  The  condition  of  the  normal  velocity  component  of  the  fluid  at  the  surface  of  the  oylin- 
dor  being  equal  to  the  normal  component  of  the  forced  velocity  of  the  cylinder,  i.e.,  if  vn  is 
the  component  of  the  forced  velocity  of  the  cylinder  in  the  direction  of  the  outgoing  unit 
normal  vector  n,  theu 


n  •  -  «  [5] 

18 

this  kinematic  boundary  condition  being  satisfied  at  the  mean  (rest)  position  of  the  cylindrical 
surface;  and 

5.  The  radiation  condition  that  the  disturbed  surface  of  the  fluid  takes  the  form  of  regular 
progressive  outgoing  gravity  waves  at  large  distances  from  the  cylinder. 

According  to  Wehausen  and  Laitone,  the  complex  potential  at  s  of  a  pulsating  point 
source  of  unit  strength  at  the  point  £  in  the  lower  half  plane  is 

<?*(«.  [log  <*  -  0  -  log  (*  -  () 

"I 

- - -  dAj  cos  tot  [0] 

-  ~  ?>  sin  to  t 

so  that  the  real  poinUsource  potential  is 

1G*  («,<;<)!  t7] 


8 


where  t  •»  a  +  iy, 

C-(  +  i>h  end 

v  m  a »a/y. 

Letting 

(?(*,/)- Be  -i  |fcg  (*- C)  -  log  <#- o 

/  -/*(*  -  C)  "I 

- —  d*J  [  8] 

“ft] 

then 

Equation  (9l  satisfies  the  radiation  condition  and  Equations  [1]  through  [4].  Another  expres¬ 
sion  satisfying  all  these  relations  is 


H 


Re  «<?(.,  C) 


[10] 


Since  the  problem  is  linear,  a  superposition  of  Equations  [9]  and  [10]  results  in 
the  velocity  potential 


•<c°  (*.y;f) 


e~i6>tde 


[11] 


where  C0  is  the  submerged  contour  of  the  cylindrical  cross  section  at  its  mean  (rest)  posi¬ 
tions  and  Q(e)  represents  the  complex  source  density  as  a  function  of  position  along  C0. 
Application  of  the  kinematic  boundary  condition  on  the  cylinder  at  s  yields 


Be  (a  •  y)  f  «(#)ff(s,C)*-0 
•&0 

let  (s  •  p)  f  Q{e)  Q{e ,  f)  de  -  A(">  u  i 

‘  Je. 


.(*) 


I 
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where  denotes  the  amplitude  of  oscillation  sad  a*"*  the  direction  cosine  of  the  normal 
velocity  at  s  on  the  cylinder.  (Both  j4(*1  and  depend  on  the  mode  of  motion  of  the  cylin¬ 
der,  as  will  be  shown  in  the  following  section.)  The  fact  that  Q($)  is  oomplex  implies  that 
Equations  [IS]  represent  a  set  of  two  coupled  integral  equations  for  the  real  functions 
Re{Q(s)\ and  /si  [<?(<)(.  The  solution  of  these  integral  equations  and  the  evaluation  of  the 
kernel  and  potential  integrals  are  described  in  the  following  section  and  in  Appendices  B 
and  C,  respectively. 


SOLUTION  OF  PROBLEM 

Take  the  *  axis  to  be  coincident  with  the  undisturbed  free  surface  of  a  conventional 
two-dimensional  Cartesian  coordinate  system.  Let  the  cross  sectional  contour  CQ  of  the 
submerged  portion  of  the  cylinder  be  in  the  lower  half  plane,  the  y  axis  being  the  axis  of 
symmetry  of  CQ.  (Since  ship  sections  are  symmetrical,  this  investigation  is  confined  to 
bodies  with  right  and  left  symmetry.)  Select  N  +  1  points  (fj,  Tjj),  (f2,  ija),  . . .  ,  (fN,  i)N),  and 
+  j,  qN  +  1)  of  CQ  to  lie  in  the  fourth  quadrant  so  that  (£,,  ijx)  is  located  on  the  negative 
y  axis.  For  partially  immersed  cylinders,  +  v  i jN+l)  is  on  the  positive  a  axis.  For  fully 
submerged  bodies,  (N  +  i  “  and  jjjv+1  <0.  Connecting  these  W  + 1  points  by  successive 
straight  lines,  N  straight  line  segments  are  obtained  which,  together  with  their  reflected 
images  in  the  third  quadrant,  yield  an  approximation  to  the  given  contour  as  nhown  in  Figure  1. 
The  coordinates,  length,  and  angle  associated  with  the  /th  segment  are  identified  by  the  sub¬ 
script  j,  whereas  the  corresponding  quantities  for  the  reflected  image  in  the  third  quadrant  are 
denoted  by  the  subscript  so  that  by  symmetry  m  -  and  ij_y  »  i\-  for  lS/S*  +  L 
Potentials  and  pressures  are  to  be  evaluated  at  the  midpoint  of  each  segment.  The 
coordinates  of  the  midpoint  of  the  ttb  segment  are 

*i  =  (fi  +  Zi  +  l)  /2»  y,-  -  (Vi  +  Vi  +  ,)  /2  [18] 

for  1  -  i-N.  The  length  of  the  tth  segment  is 

l»,-l  “  >/(£,•  + 1  -  £<)a  +  (ii  + 1  -  n,)2  [1*1 

while  the  angle  made  by  the  tth  segment  with  the  positive  x  axis  is  given  by 


■  tan~ 1 


+  1  “  ’I  i 


[15] 
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The  outgoing  unit  vector  normal  to  the  cross  section  at  the  tth  midpoint  (*t,  yj)  is 


-  i  sin  |  cos  mi 


U«J 


where  i  end  |  ere  twit  vectors  in  the  directions  of  increasing  s  and  jr,  respectively. 

The  cylinder  is  forced  into  simple  harmonic  motion  with  radian  frequency  at,  according 
to  the  displacement  equation 


g<*0 •  eoe  01 1 


[171 


for  m  ■  8,  I,  or  4,  corresponding  to  sway,  heave,  or  roll,  respectively.  The  rolling  oscillations 
are  about  an  axis  through  a  point  (0,  is  the  symmetry  plane  of  the  cylinder. 

In  the  translational  modss,  any  point  of  ths  cylinder  moves  with  velocity 


for  sway,  or 


■  -  i  iO)  •  sin  at 


*>  sin  at 


(181 


[191 


for  heaving  oscillations. 

Rolling  motion  about  (0,  y0)  la  illustrated  in  Figure  8.  Considering  s  point  y$  > 
co>  an  inspection  of  Figure  8  yields 


-  A2  +  (»<  -  vj* 


8,  -  tan  1 -  a  am  1  — — —  ■  cos  1  - 

‘  #,  Bt  B, 


Therefore,  by  elementary  two-dimensional  kinematics,  the  unit  vector  in  the  direction  of  in¬ 
creasing  8  is 


t{  m  -  i  sin  6t  +  {  cos  9t  ■  -  ■ 


*<-*0  *< 


- •*»»**.. 


so  that 


v<4>  -  B.  <S<4>  m  O)  d<4)  [(y4  -  y0)  I  -  wi  I  ]  sin  at 


m 


The  normal  components  of  the  velocity  e/M)  -  »4-  •  v^"}  tt  the  midpoint  of  the  i th 
segment  (»4,  y{)  ere 


•  -  e*  sin  a(  sin  rnt 
*/3>  »  a  cos  m4  sin  at 
»/4)  -  ai  4<4>  [(y4  -  y0)  sin  a{  +  *4  cos  «4]  sin  at 
far  the  swaying,  heaving,  and  rolling  modes,  respectively.  Defining 


[21] 


nfm) ' 


a,  sin  cut 


then,  consistent  with  the  previously  mentioned  notation,  the  direction  cosines  for  the  three 
modes  of  motion  are 

n/2)  -  -  sin  a  ( 


»/3)  -  cos  m4  [22]  j 

I 

«/4>  -  (yf  - y0)  8in  •*  +  *,• 008  j 

i 

i 

Equations  [2£]  illustrate  that  heaving  is  symmetrical,  i.e.,  nj3.>  »  n/3).  Swaying  and  rolling, 
however,  are  antisymmatrical  modes,  i.e.,  -  -  nfa\  and  ni4)  -  -  nf4\ 

Equations  [12]  are  applied  at  the  midpoints  of  each  of  the  N  segments,  and  it  is 
assumed  that  over  an  individual  segment  the  complex  source  strength  Q(a)  remains  constant, 
although  it  varies  from  segment  to  segment.  With  those  stipulations,  the  set  of  coupled 
integral  Equations  [12]  becomes  a  set  of  2N  linear  algebraic  equations  in  the  unknowns 

Re  ie<*>  (*;)!  -  Qfm>  and  Im  !$<«>  {»f)  |  -  Qj*) 


7 


[28] 


Thus,  for  i  •»  1,  2,  . . .  ,  N 


N  N 

/xi  <?/->//;),  yxt  QtiViW- o 


/V 


yv 


-  x  c/">4">  + 


/-i 


,Jt  «nVi/ 


/<*>  -  w  4<w)  »<"> 


where  the  superscript  (■>)  denotes  the  mode  of  motion.  The  “influence  coefficients*' 
end  <f^M>  and  the  potential  yf,  l)  are  evaluated  in  the  appendix.  The  resulting 

velocity  potential  consists  of  a  term  in  phase  with  the  displacement*  and  a  term  in  phase 
with  the  velocity. 

The  hydrodynamic  pressure  at  along  the  cylinder  is  obtained  from  the  velocity 

potential  by  means  of  the  linearised  Bernoulli  equation 

y<i,)  (*/. Vi, *)  *  - p  *em> (»(,  *»;  0 


as 


/m)  (*p  y0  a;  t)  -  p<m)  (»,,  yt\  a)  cos  at  +  p  <"  >  y,;  a)  sin  at  [25] 

where  pja)  and  are  the  hydrodynamic  pressures  in  phase  with  the  displacement  and  in 
phase  with  the  velocity,  respectively,  (p  denotes  the  density  of  the  fluid  in  Equation  [24]). 

As  indicated  by  the  notation  of  Equations  [241  and  [25],  the  pressure  as  well  as  the  potential 
is  a  function  of  the  oscillation  frequency. 

The  hydrodynamic  force  or  moment  (when  •»  -  4)  per  unit  length  along  the  cylinder 
necessary  to  sustain  the  oscillations  is  the  integral  of  pim)  ♦  n(s>>  over  the  submerged  contour 
of  the  cross  section  CQ.  It  is  assumed  that  the  pressure  at  the  ith  midpoint  is  the  mean  pres¬ 
sure  for  the  ith  segment,  so  that  integration  reduces  to  summation,  whence 

M 

Jf(*°  («)  -  2  2^  pj",)  a)  »/") 

and  ' 

*<->(„). 2  pr>(«i.y<;«)»/",> 

for  the  added  mass  and  the  damping  forces  or  moments,  respectively. 


•Stott  author*  rater  to  this  Ufa  at  a  coaapooant  la  pkata  with  tht  account!  oo,  Howarar,  da*  to  dlapUcwnant 
Equation  [l7],  thla  author  Saaan  it  eor*  appropriate  to  rater  to  thia  tana  as  halo*  ISO  de*ree«  out-of-phase  with 
tho  accalaratioo  0*  Ia-phaso  with  tho  dloplacaawat. 


[26] 

[27] 
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The  velocity  potentials  for  very  smell  end  very  large  frequencies  ere  derived  end 
discussed  ia  Appendix  D. 

LIMITATIONS  OF  GREEN'S  FUNCTION  -  INTEGRAL  EQUATION  METHOD 

F.  John13  proved  the  existence  ead  uniqueness  of  the  solutions  to  the  three*  sad  two- 
dimensional  potential  problems  pertaining  to  oscillations  of  rigid  bodies  ia  e  free  surface. 
The  solutions  were  subject  to  the  provisions  that  no  point  of  the  immersed  surface  of  the 
body  would  be  outside  a  cylinder  drawn  vertically  downward  from  the  intersection  of  the  body 
with  the  free  surface  and  that  the  free  surface  would  be  intersected  orthogonally  by  the  body 
in  its  mean  or  rest  position. 

John  also  showed  that  for  a  set  of  discrete  “irregular”  frequencies  the  Green’s 
function-integral  equation  method  failed  to  give  a  solution.  He  demonstrated  that  the 
irregular  frequencies  occurred  when  the  following  adjoint  interior-potential  problem  had 
eigenfrequencies.  Let  p(*,y)  be  such  that 

1.  <lixx  +  $yy  »  0  inside  the  cylinder  in  the  region  bounded  fay  the  immersed  surface  of 
the  body  and  the  extension  of  the  free  surface  inside  the  cylinder; 

2.  -  vk  i/f  «  0  on  the  extension  of  the  free  surface  inside  the  cylinder,  vk  being  the 
wave  number  corresponding  to  the  irregular  frequency  k  -  1,  2,  8, ...  ; 

8.  ^  -  0  on  the  surface  of  the  cylinder  below  the  free  surface. 


For  a  rectangular  cylinder  with  beam  B  and  draft  T  (see  Figure  8),  the  irregular  wave 
numbers  may  be  easily  obtained  by  separation  of  variables  in  the  Laplace  equation.  Sepa¬ 
rating  variables  gives  the  eigenfunctions 

i/ik  m  Bk  sin  (k  w  w/B)  sinh  (A  w  T/B) 

for  k  «  1,  2,  8, ...  ,  etc.,  where  the  are  Fourier  coefficients  to  be  determined  from  an 
appropriate  boundary  condition.  Applying  the  free-surface  condition  (2)  on  y  -  T  for 
0  <  a  <  B,  the  eigenwave  numbers  (or  irregular  wave  numbers) 

vk  «  (k  w/B)  ooth  (k  w  T/B)  (28] 

are  obtained  for  k  -  1,  2,  ... ,  etc.  In  particular,  the  lowest  such  irregular  wave  number  is 
given  by 


-  (w/B)  coth  (*  T/B) 


[29] 
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Keeping  T  fixed  in  Equation  [89]  but  letting  B  vary  and  setting  4  -  w/B,  then  from  the 
Taylor  expansion 

b  ooth  (»f)  -  4  [l/(»r>  +  bT/9  -  (5T)3/ 48  + . . .] 

it  is  seen  that  as  b  ■*  0,  which  is  equivalent  to  B  +  •»,  vx  -  1/f.  Therefore,  for 
rectangular  cylinders  of  draft  T, 


vx  *  1/f  [80] 

a  relation  that  John  proved  for  general  shapes  complying  with  the  restrictions  previously 
outlined.  Figure  4  depicts  the  first  irregular  wave  number  for  rectangular  cylinders  as  a 
function  of  the  beam-to-draft  ratio,  keeping  the  draft  constant.  For  a  beam-to-draft 
ratio  of  8.5,  vx  «  1.48;  while  for  B/T  -  2,  vx  -  1.71. 

At  an  irregular  frequency  the  matrix  of  influence  coefficients  of  Equations  [88]  be* 
comes  singular  as  the  number  of  defining  points  per  cross  section  increases  without  limit, 
i.e.,  as  N  *  •».  In  practice,  with  finite  N,  the  determinant  of  this  matrix  becomes  very  small 
not  only  at  the  irregular  frequency  but  also  in  an  interval  about  this  frequency.  This  interval 
caa  be  reduced  by  increasing  the  number  of  defining  points  JV  for  the  cross  section.  Figures 
6,  8,  and  7  illustrate  this  phenomenon  and  its  effect  on  the  calculated  added  mass  and  damp* 
ing  coefficients. 

These  results  are  significant  in  the  application  to  ship  sections.  Most  surface  vessels 
have  nearly  constant  draft  over  the  length  of  the  ship,  and  the  maximum  beam  occurs  at  or 
near  midship,  where  the  cross  section  is  usually  almost  rectangular,  so  that  for  most  surface 
ships  the  first  irregular  frequency  <u,  is  less  for  the  midsec ti on  than  for  any  other  cross 
section.  For  a  ship  with  a  7:1  length-to-beam  and  a  5:2  beam-to-draft,  the  first  irregular 
wave  encounter  frequency-in  nondimeosional  form  with  L  denoting  the  ship  length-occurs  at 

e>x  'Jh/g  *  5.09 

which  is  beyond  the  range  of  practical  interest  for  ship-motion  analysis.  Therefore,  for 
slender  surface  vessels,  the  phenomenon  of  the  first  irregular  frequency  of  wave  encounter 
is  not  too  important. 
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NUMERICAL  RESULTS  AND  CONCLUSIONS 


A  computer  program,  based  on  the  analysis  of  this  report,  was  developed,  and  calcula¬ 
tions  for  cylinders  of  various  cross  sections  were  performed  by  the  IBM  7000  electronic 
digital  computer  at  the  Center.  Two  options  in  the  program  provide  for  pressure  and  added 
mass  or  moment  calculations  for  the  asymptotic  cases  of  very  small  and  very  large 
frequencies. 

The  most  time  consuming  computer  operations  occur  when  evaluating  the  principal 
value  integrals.  These  integrals  are  of  the  complex  exponential  integral  type  and  are  cal¬ 
culated  as  power  series  in  vr,  with  r2  =  (*,  -  fy)2  +  (y(  +  q;)2,  and  v  »■  u2/g.  These  series 
(derived  in  the  appendix)  converge  fairly  rapidly  for  partially  submerged  shiplike  sections  at 
frequencies  lower  than  the  firBt  irregular  frequency.  This  results  in  an  average  computer 
time  of  20  seconds  per  frequency  for  half  sections  defined  by  21  input  points.  However,  for 
fully  immersed  bodies,  r  is  considerably  larger,  so  that  the  much  slower  convergence  of  the 
series  increases  the  computation  time  greatly.  For  a  given  contour  and  frequency  the  com¬ 
puter  time  varies  as  the  square  of  the  number  of  input  points,  which  is  limited  to  a  mesh 
fineness  of  46  points  per  half  section. 

The  frequency  input  to  this  program  is  in  the  form  of  wave  numbers,  nondimensional- 
ized  by  some  principal  linear  dimension  such  as  the  draft  or  the  half  beam.  The  output  of 
this  computer  program  includes  the  pressures  in  phase  with  the  displacement  and  in  phase 
with  the  velocity  at  each  of  the  segmental  midpoints,  the  added  mass  or  added  moment  and 
damping  coefficients  in  normalized  form,  and  a  scaled  version  of  the  determinant  of  the 
matrix  of  influence  coefficients. 

For  partially  submerged  cylinders  the  added  mass  or  moment  and  the  damping  are 
scaled  by  &>2  times  the  amplitude  of  oscillation  times  the  displaced  mass  or  displaced 
moment  of  inertia  of  a  semicircular  cylinder  of  radius  equal  to  the  principal  dimension.  For 
fully  immersed  bodies  these  quantities  are  normalized  by  u2  times  the  oscillation  amplitude 
times  the  displaced  mass  or  displaced  moment  of  inertia  of  a  circular  cylinder  of  radius  equal 
to  the  principal  dimension.  It  is  assumed  in  all  cases  that  the  cylinders  are  oscillating  with 
unit  amplitude 

The  graphs  of  Figures  5  through  26  were  obtained  with  the  aid  of  the  SC  4020 
Charactron  plotter,  an  important  piece  of  peripheral  equipment  to  the  IBM  7090  digital 
computer.  These  figures  illustrate  the  behavior  of  the  hydrodynamic  pressure  and  normalized 
force  components  as  functions  of  the  nondimensional  wave  number  for  various  cylindrical 
shapes.  In  Figures  8  through  18,  the  added  mass  and  damping  coefficients  as  well  as  the 
hydrodynamic  pressure  components  are  plotted  against  <u2  B/(2g)  for  the  heaving  and  swaying 
semi -immersed  circular  cylinder,  defined  by  21  equally  spaced  input  points  around  the 
half  section.  The  results  for  the  heaving  mode  agree  well  with  the  Porter  results.  Figure  14 
depicts  the  added  mass  and  damping  coefficients  of  a  fully  submerged  heaving  circular  cylin¬ 
der,  defined  by  21  equally  spaced  points,  whose  axis  is  submerged  to  a  depth  of  1.25  times 
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its  radius.  The  values  obtained  in  this  case  are  in  good  agreement  with  the  results  obtained 
by  Ogilvie  from  his  first  order  theory.  Figures  IS  through  17  show  the  added  mass,  damping, 
and  pressures  of  a  partially  immersed  heaving  ogival  cylinder  with  draft/diameter  equal  to 
9/10  and  beam/diameter  equal  to  8/5  as  functions  of  the  nondimensional  wave  number 
<u2  B/(2g).  This  is  an  instance  of  a  bulbuous  section  similar  to  the  cylinders  discussed  by 
Motors  and  Koyama, 14  which  have  the  property  that  the  damping  vanishes  at  some  frequency. 
Figure  18  indicates  aero  damping  at  about  u2B/{ 2p)  *  0.7. 

Figures  18  through  20  illustrate  the  added  mass,  added  moment,  and  damping  coeffi¬ 
cients  for  all  three  modes  of  oscillation  of  a  cylinder  whose  cross  section  is  in  the  shape  of 
the  midsection  of  a  Series  60,  Block  70  ship  defined  by  16  input  points  per  half  section  as 
well  as  the  same  parameters  for  a  geometry  defined  by  only  7  input  points.  The  close  agree¬ 
ment  between  these  two  sets  of  output  data  (within  10  percent  of  each  other)  suggest  the 
possibility  of  using  a  coarse  sectional  input  mesh  for  the  application  of  this  method  to  the 
strip  theory  of  ship  motions.  The  computing  time  for  the  16  point  geometry  was  roughly  six 
times  that  for  the  7  point  half  section.  The  graphs  in  Figure  18  are  very  similar  to  the  curves 
for  the  added  mass  and  damping  coefficients  of  the  heaving,  full-form,  Lewis-type  cylinder 
of  Porter.  Figures  21  through  26  show  the  Bway,  heave,  roll,  added  mass  or  moment  and 
damping  coefficients  for  the  bow  and  stern  sections  of  the  same  ship,  defined  fay  7  points  per 
half  section.  Observe  that  the  roll  damping  for  these  sections  is  considerably  larger  than  for 
the  midsection,  being  largest  for  the  extremely  slim  bow  section. 

Recall  that  three  assumptions  were  made  to  obtain  numerical  solutions  to  the  boundary- 
value  problem.  First,  it  was  assumed  that  the  cylindrical  cross  section  could  be  approximated 
by  a  polygonal  succession  of  straight-line  segments  defined  by  a  finite  number  of  points. 
Second,  the  hypothesis  of  the  constancy  of  the  source  strengths  along  an  individual  segment 
was  employed  to  convert  the  integral  equations  into  a  set  of  linear  algebraic  equations. 

Third,  the  supposition  of  mean  hydrodynamic  pressures  over  the  individual  segments  made  it 
possible  for  the  integrations  of  these  pressures  to  be  transformed  into  summations.  However, 
as  the  number  of  segments  is  made  to  increase  indefinitely,  while  the  largest  of  these  seg¬ 
ments  converges  to  zero  length,  the  polygonal  approximation  approaches  the  cross  sectional 
contour,  and  the  source  strengths  and  pressures  become  continuous  functions  of  position  (and 
frequency)  along  this  contour. 

The  Green's  function-integral  equation  method  os  described  in  this  paper  is  applicable 
to  any  two-dimensional  simply  connected  shape  (even  those  with  severe  corners,  bilge  keels, 
bulbs,  fins,  etc.),  whereas  the  mapping  of  those  odd-shaped  contours  onto  the  unit  circle 
requires  a  great  number  of  terms  in  the  Thecdorsen  transformation.  Smith  reported  the 
necessity  of  obtaining  16  Theodorsen  coefficients  for  the  mapping  of  a  bulbous  form  with  an 
expenditure  of  8  minutes  of  computer  time. 
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APPENDIX  A 


EVALUATION  OF  PRINCIPAL  VALUE  INTEGRAL 


The  real  and  imaginary  parts  of  the  principal  value  integral 


*°*  e~ik  & 


are  used  in  evaluating  some  of  the  kernel  and  potential  integrals. 

The  residue  of  the  integrand  at  It  ■  vis  «“*  **(*  •  (),  ao  that 


e-ik  (*  -O 


g~  (*  -  i) 


dh^ine-iv^  "  £>  [31] 


where  the  path  of  integration  is  the  positive  real  axis  indented  into  the  upper  half  plane  about 
k  ■  v. 

Note  that  v  -  <ua/g  >  0,  !m  »  <  0,  and  /«  <  <  0.  The  transformation  «-£(*-  v)  (g  -  £) 
converts  the  contour  integral  on  the  right  side  of  Equation  [81]  to 


- - - -  I  — 

J  V-*  I  -  U> 

o  *-  <v  <i  -  o 


e-‘v(i  -  o  £-i  [-^(a  -  ft] 

-  «“  iv<*  “  ^  jy  +  log  Mv(s  -  01 

+  Jp  (-l)n[-*v(»-C)H 

n  »  1  ' 


where  y  •  0.5772  ...  is  the  well  known  Euler-Maacheroni  constant.  (See  Abramovitz  and 
Stegun15  for  the  definition  of  Ev) 

Setting  r  -  |-tV(a  -  ft  I  and 

8  -  tan-1  [f*»(-iv(a  -  ft)/ffe(-tv(a  -  ft)]  +  ir 
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the  following  expression  is  obtained  for  Equation  [81] 

/**  *-**  <» «  Z) 


•in  v(*  -  £)] 


*il  v  + 


— 

1 °*'Z 


r"  cos  (n  0) ' 


Separating  Equation  [03]  into  its  real  and  imaginary  parts  yields 


M 

»  V  f  «<r  +  1,5  *  oos  *(*  - 

J  - «>^M)[C(r,D 


oos  v(»  -  f) 


*  S(r,  fl)  sin  v(»  -  £)] 


pv  f  «<),  +  ,»*sin*(«-f) 

J  - TT*-* - -  + 


•in  v(a  -  () 


~  S(r,  9)  oos  v(»  -  £)] 


provided  that 


,«)-y  +  10g  r  +  yjl'Q*  < 
m  n  •  n 


r"  oos  (n  3) 


S(r,  9)  -  g  +  _  ,in  <"3 

^  a  •  nl 
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sinco 


(+iri 

d  C  ■  d(  +  idt) 

*  da  coa  * .  +  id*  sin 
-timid* 

for  (  along  the  ;'th  segment.  Therefore,  da  »  e  d(,  and 


R*  !(•»/  •  9)  f  log  (*  -  0  da  l 

l"  **  */  *  l#/  log  <*  ~  j  I 

i "" 


p  tX.  1o«  <a  -  di 

Setting  C  «  s  -  £,  the  last  integral  becomes 


a 

i  r~c> 

Re 

J 

*  logr^C' 

O 

[87] 


-  wg  (*/  -  C/)  -  *rg  (*/  -  C/+  j  )  -  » 

If  i  d  j,  differentiation  under  the  integral  sign  may  be  performed,  so  that 


|(»«  •  ?)  j  log  (s  -  O  da  J  j 

Vr  (*«-  f;)2  +(y<-’?/)2 

- - - - 

(*<-  ^/+i) 

f  .  y»  -  *7/  Vi-^z+n 

+  cos  («.  -er,)  I  tan-*  - —  -  tan  1  - - — 

4  '  L  *<-£, 


[881 
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For  the  integral  containing  the  log  (a  -  2)  term,  da  •  e  *  d%,  so  that 

Be  jflog  (a-?)*#  |  } 


*«* », 


•  sin  (< 


Vr  +(y,  +  i/)i 


'(*i-l/+1)2  +  (y<  +  »i/+i)a 

r  _i  y<+,»/  _t  y<+,»/+n 

♦  [u.-  ^-u.-  ^7,J 


[Ml 


f  f  ? 

*«  J<"<F>  J  *P.F.  | 

- __  dfc  l 

v-*  •-*<{ 

1r  .  -<*<*<-  o 

+-/)J_  di  ^  p.f.J^  rry 

j  -W 


t40] 


-p.r. 


-  cos  («4  +« 


■>  [""J 


S^^sin*^- 


€j) 


-  P.K. 


I 


v  -  * 

e‘(r^u/+i)8i 


rf* 


sin  *(*<  -  £;>1) 


A 


*] 


he  first  integrsl  on  the  right  side  of  'Equation  [86]  becomes 

k  -F)  J  e-^'-i'dal  \ 

l  *t  !*-*</ 

>  -  sin  (*(  +  a  j)  cos  -  (,)  -  /(y<+  ^  cos  ,(*<-f/+l)] 

■  cos  (<■<+«,)  ^l/(y<  +  ,/>8in  i/(*<-f/)-«^r<  +  ,/  +  1^sin 


[41] 


l 


The  kernel  integrals  over  the  image  segments  are  obtained  from  Equations  138]  through  [41] 
by  replacing  £/  +  1,  end  tt.  with  £_;  =  -  ij,  £_  (/  + -  -  f/  +  i>  and  a  -/  “  _  “/’ 
respectively. 


APPENDIX  C 


EVALUATION  OF  POTENTIAL  INTEGRALS 


The  velocity  potential  for  the  aith  mode  of  oscillation  at  the  ith  midpoint  (ari ,  y()  is 


i  £  «,*«  { J  [' 


log  (>i  -<)-  log  (a,  -  C) 


* 2 ^ v  J  -TT*—"!* 

0  J 

-  (-  1)*  J  j^log  (a,  +  ?)  -  log  (a,  +  O 


+2PK-  J  d‘\ 


/•  1  Way 


cos  <uf 


ala 


-Mr  j  .-'-‘*<*0*1 


•  sin  <ut 


(421 


-/ 


The  integration  of  the  log  (i.  -  £)  term  is  straightforward,  yielding 


"(/  log  (*<  -  O  <*«J 


cos  a 


7  £<#,  -  (j)  log  \/(*j  -  (,)%  +  (y,  -  Vj)3  +  £,  -  £,  +  i 
« (*<  -  €j  +  j)  log  >/(*,•  -£/+l)2 +(y<-'i/+1)2 

,  (Ti-»/)  ,  (Vi- »»/+1n 


[431 
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In  the  integration  of  the  lop  (*  -  Z)  term,  note  that  ij?-  and  ij^  +  1  are  replaced  by  -  t/.  and 
-  jj y  +  1>  respectively. 

To  evaluate  the  potential  integral  containing  the  principal  value  integral,  proceed  in 
the  following  manner.  For  an  arbitrary  a  in  the  fluid  domain 


where  the  change  of  integration  is  permissible  since  only  one  integral  requires  a  principal 
value  interpretation.  After  dividing  by  v  and  multiplying  by  u  -  k  +  k  under  the  integral 
sign,  the  last  expression  becomes 
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Regarding  the  first  integral  in  Equation  [44]  as  a  function  of  a, 


0 


[45] 


Differentiating  Equation  [45]  with  respect  to  a  gives 


F  '(a)--f 


-ik  (* 
e 


00 

J 

0 


-Ik  <*  - 
e  ' 


1  1 

a-2,  +  1 


so 


F(a)  -  log  (a  -  <y)  -  log  (a  -  ?/+  ,)  +  K  [46] 

where  K  is  a  constant  of  integration  to  be  determined  presently.  Since  F  (a)  is  defined  and 
analytic  tor  all  a  in  the  lower  half  plane,  and  since  by  Equation  [45],  lim  ^  ^  F (a)  -  0,  it 
follows  from  Equation  [46]  that  K  -  0.  Therefore, 


Re 


mRe 


“  -ik  (*i  -  £) 
daP.V.  J  - - — - 
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APPENDIX  D 


VELOCITY  POTENTIALS  FOR  VERY  SMALL  AND  VERY  LARGE  FREQUENCIES 


For  very  smell  frequencies,  i.e.,  es  <u  ■*  0,  the  free -surface  condition  Equetion  [8} 
of  the  section  formulating  the  problem  degenerates  into  the  wall-boundary  condition 

♦y -0  [49] 

on  the  surface  of  the  fluid  outside  the  cylinder,  whereas  for  extremely  large  frequencies, 
i.e.,  when  <u  •»,  the  free-sorface  condition  becomes  the  “impulsive”  surface  condition 

«  «  0  [50] 


ob  |r  ■»  0  outside  the  cylinder. 

Equations  [S3,  [4],  and  [51  remain  valid  for  both  asymptotic  cases.  The  radiation 
condition  is  replaced  by  a  condition  of  boundedness  at  infinity. 

Therefore,  there  is  a  Neumann  problem  for  the  case  <■>  •»  0  and  a  mixed  problem  when 
The  appropriate  complex  potentials  for  a  source  of  unit  strength  at  a  point  C  in  the 
lower  half  plane  are 


o  -  D°«  [•-*>♦  1°S  (•  -  ?>3  +  *0  foH 

and 

C-(*’ [log <)]  +  *..  [#21 

for  the  Neumann  and  mixed  problems,  respectively,  where  K0  and  Km  are  constants  not  yet 
specified. 

Let 

(*»  r,  l  n)  -  [o,  (*,  O) 

so  that  the  velocity  potentials  for  the  m  th  mode  of  motion  are 

•i")  (*.  y)  -  J  Q,<")  (•)  4>m  (*.  y;  & »»)  Ml 

c* 


for  a  -  0,  and  a  •  <»,  where  is  the  expression  for  the  source  strength  as  a  function  of 
position  along  the  submerged  contour  of  the  cross  section  CQ.  An  analysis  similar  tc  the  one 
in  the  section  on  formulating  the  problem  leads  to  the  integral  equation 

("  •  r)  J  Q^m)  (»)  4>a  (*,  y;  £.  q)  *  -  »{m)  [Ml 

co 

which,  after  application  at  the  N  segmental  midpoint,  yields  a  set  ot  N  linear  algebraic 
equations  in  the  N  unknown  source  strengths  Qj. 

It  remains  to  be  shown  whether  these  two  problems  are,  in  the  language  of  potential 
theory,  well  posed,  i.e.,  whether  the  solutions  to  these  problems  lead  to  unique  forces  or 
moments.  The  mixed  problem  raises  no  difficulty,  since  a 8  »■*•»,  Om  (s,  C)  *  0.  In  fact 
Afw  -  0,  which  can  be  inferred  from  the  pulsating  source-potential  Equation  [8]  by  letting 

V  +  O0, 

Considering  the  Neumann  problem,  note  that  the  constant  KQ  in  the  Green’s  function 
Equation  [51]  yields  by  integration  an  additive  constant  AT  to  the  potential.  However,  for  a 
completely  submerged  cylinder  the  cross  sectional  contour  C0  is  a  simply  closed  curve,  so 
that  the  contribution  of  K  in  integrating  the  product  of  the  pressure  with  the  direction  cosine 
of  the  body-surface  velocity  vanishes.  For  partially  submerged  bodies  C0  is  no  longer 
closed.  But  since  nf"*  •  -  for  m  being  even, 

J  Afn<">ds«  0 

co 

so  that  the  swaying  force  and  rolling  moment  are  unique.  The  heaving  foroe  on  a  partially 
submerged  cylinder  is  not  unique  for,  in  this  case,  «  n/3),  no  that 

e 

J  K  dt  to 

Co 

The  constant  Afc  may  be  obtained  fay  letting  v  *  0  in  the  pulsating  source-potential 
Equation  [8]. 
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|(fl'9l) 

Fignra  1  -  Polygonal  Approximation  to  Immersed  Part 
of  Cylindrical  Cross  Section  Cn 


(21  Input  Points)  (Defined  by  21  Input  Points  pet  Half  Seettou) 
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(7  hprt  Petals) 


ire  Added  Maas  end  Damping  Coefficients  Figure  25  —  Sway  Added  Mass  and  Damping  Coefficients  for 

action  of  Series  60,  Block  70  Ship  for  Aft  Section  of  Series  60,  Block  70  Ship 


Figure  96  -  Roll  Addad  Moment  and  Damping  Coefficient*  for  Aft  Section 
of  Seriee  60,  Block  TO  Ship 
(7  hfit  Palate) 
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